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INTRODUCTION 


'^Projectile  penetration 


studied* hy  military  engi- 


neers since  ancient  times,  but  the  analytical  inve^igation  of  this 
phenomenon  was  not  undertaken  until  the  eight eenth^entury.  Robins 
(1),  Euler  (2),  and  Poncelet  (3)  developed  the  firs^|mpirical 
equations  describing  the  axial  force  which  opposes  the  motion  of  a 
rigid  projectile  as  it  passes  through  a.  deformable  solid  target^i;**^^ 
Experimental  data  consisting  of  iinpact  velocity  and  final  depth  of 
penetration  were  used  to  estimate  the  coefficibnts  in  the  relation 
between  axial  force  and  projectile  velocity.  However,  these  coef- 
ficients were  not  specified  in  terms  of  standard  material  proper- 
ties, and  so  they  were  applicable  only  for  the  particular  targets 
for  which  they  were  determined. — 

/ 

Since  the  mid-twentieth  century,  the  availability  of  high- 
speed digital  computers  has  made  it  possible  to  analyze  the 
projectile/target  interaction  by  means  of  two-dimensional  finite- 
difference  techniques.  The  constitutive  relations  for  both  the  tar- 
get and  the  projectile  can  be  accurately  specified,  and  many  of  the 
details  in  the  penetration  process  can  be  realistically  simiilated. 
Analyses  of  this  kind  are  valuable  because  of  the  detailed  informa- 
tion which  they  provide,  but  they  are  usually  too  costly  and  cumber- 
some for  large  scale  parameter  studies. 

In  recent  years,  significant  progress  has  been  made  toward 
the  deve].opment  of  a simple  relation  between  penetration  resistance 
and  fundamental  target  properties  such  as  density  and  compressive 
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strength.  Bishop,  Hill,  and  Mott  (U)  analyzed  the  quasistatic  ex- 
pansion of  spherical  and  cylindrical  cavities  in  metals  and  observed 
a correlation  between  the  resistance  to  cavity  expansion  and  the 
resistance  to  indentation  by  a rigid  static  punch.  Goodier  (5)  used 
an  analogy  with  the  dynamic  'eXpaHsion  of  sphericeLL  cavities  in 
metals  to  develop  a penetration  theory  for  rigid  spherical  projec- 
— tiles .' ^In  the  present  work,  the  quasistatic  expansion  of  a spheri- 
cal cavity  in  a concentrically  layered  medium  will  be  analyzed  in 
order  to  obtain  an  approximation  for  the  radial  stress  at  the  sur- 
face of  a slowly  expanding  spherical  cavity  in  a vertically  layered 
medium  with  a plane  interface.  This  approximation  will  then  be  used 
to  construct  a simple  equation  of  motion’  for  a rigid  cylindrical 
projectile  (with  a conical  or  ogival  nose  shape)  which  penetrates  a 
layered  target. 


SPHERICAL  CAVITY  EXPANSION  IN 
A CONCENTRICALLY  LAYERED  MEDIUM 


Consider  a slowly  expanding  spherical  cavity  surrounded  by 
an  infinite  medium  composed  of  two  distinct  concentric  layers 
(Figure  l).  The  materials  surrounding  the  cavity  exhibit  elastic- 

plastic  behavior  without  vol- 
ume change  (Figure  2).  The 
geometry  is  spherically  sym- 
metric and  the  material  equa- 
tion of  equilibrium  is 


(1) 


where  r is  the  radial  coor- 
dinate and  Oj.  and  oq  are 
respectively  the  radial  and 
circumferential  components  of 
stress.  (Compressive  stress 
and  strain  are  taken  to  be 
positive. ) The  stress-strain 
relation  for  the  material  in 
the  elastic  state  is 
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Figure  1.  Slowly  expanding 
spherical  cavity  in  a concen- 
trically layered  medium 


"r  - "e  = 3 


and  in  the  plastic  state  it  is 
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Figure  2,  Material  response  to  uniaxial  stress 


“r  - ”8  “ * 


(3) 


where  E is  the  elastic  modulus,  Y is  the  compressive  strength, 
and  Sj.  and  £9  are  respectively  the  radial  and  circumferential 
strains.  Both  materials  are  incompressible,  and  the  relation  be- 
tween and  £9  is 


^ +229  = 0 


ih) 


Equations  2 and  4 are  now  combined  so  that  the  stress-strain  rela- 
tion in  the  elastic  state  becomes 

ar  - Og  = -2E£g  (5) 

The  circumferential  strain  at  any  radial  position  is  given  by 

o 

where  ro  is  the  initial  radial  position  of  a material  particle 
which  has  been  displaced  to  position  r , Assuming  that  the  spheri- 
cal cavity  expands  from  a small  initial  radius  a^  to  an  instantan- 
eous radius  a >>.ag  , the  principle  of  conservation  of  mass  re- 
quires that 


3 3 3 3 ^3 

r - r = a - a a 
o o 

When  Equation  T is  rearranged,  it  follows  that 


(7) 
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and  the  expression  for^ the  circuffli^rential  strain  hecomes 

(‘■S 


(8) 


e0«|ln 


(9) 


Restricting  attention  to  small  strains  in  the  elastic  state,  Equa- 
tion 9 is  further  approximated  hy 


e 

0 3 ^3 

and  the  stress-strain  relation  in  the  elastic  state  now  "becomes 


(10) 


a - 
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Combining  Equation  1 alternately  with  Equations  3 and  11,  the  mate- 
rial equation  of  equilibrium  in  the  plastic  state  is 


dr 

and  in  the  elastic  state  it  is 


dcr  • }.  3 

dr  3 


(12) 


(13) 


Each  material  layer  may  contain  distinct  elastic  and 
plastic  regions  (Figure  3)  which  are  separated  by  spherical  plastic 
"fronts”  located  at  r = b and  r = k in  layers  1 and  2,  respec- 
tively. It  is  obvious  from  Figure  3 that  b £ h since  the  plastic 
region  in  layer  1 cannot  extend  beyond  the  location  of  the  layer 
interface  at  r = h . Similarly,  it  follows  that  k ^ h in  layer  2. 
Evaluating  Equations  12  and  13  at  r = b and  r = k and  equating 
the  right-hand  sides  at  each  of  these  locations,  then  it  is  found 
that 


b = a 


’2E^ 

,3^ 
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expanding  spherical,  cavity  is  centered  on  the  z-axis  at  a distance 
IIq  from  the  left  side  of  the  layer  interface.  A first-order  approx- 
imation for  the  radial  stress  a on  the  right  hemisphere  of  the 
cavity  surface  can  be  obtained  by  regarding  the  vertical  layers  as 
quasi-concentric  layers'-'vith  the 'initial  position  of  the  layer  inter- 
face located  at  a radial  distance  r = ho  from  the  center  of  the 
cavity.  The  effective  radial  position  of  the  layer  interface  is 
then  given  by  Equation  l8,  and  the  first-order  approximation  for  o 
is  given  by  Equation  20. 

PROJECTILE  EQUATION  OF  MOTION 

IN  A VERTICALLY  LAYERED  MEDIUM 

Projectile  penetration  and  spherical  cavity  expansion  rep- 
resent geometrically  dissimilar  processes.  Nevertheless  it  has  been 
observed  that  the  axial  force  which  opposes  penetration  is  more  or 
less  proportional  to  o at  low  velocities  in  materials  which  are 
composed  of  a single  semi-infinite  layer  (4,5).  It  has  been  further 
observed  that  the  relation  between  final  depth  and  impact  velocity 
is  approximately  linear  for  deep  penetration  (5),  which  implies  a 
possible  linear  relation  between  the  axial  resisting  force  and  the 
projectile  velocity. 

These  observations  will 
now  be  used  as  guidelines 
in  the  formulation  of  an 
ad  hoc  expression  for  the 
effective  stress  nomal 
to  the  frontal  surface  of 
a penetrating  projectile. 


Consider  a 

rigid  axisymmetric  pro- 
jectile which  penetrates 
a semi-infinite  target 
composed  of  two  distinct 
vertical  layers  (Fig- 
ures 5 and  6).  The  pro- 
jectile axis  of  symmetry 
coincides  with  the  direc- 
tion of  motion,  which  is 
normal  to  the  target  sur- 
face. The  quantity  a 
contains  the  pertinent 
dimensions  and  material 
properties  of  the  target 
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Figure  5.  Projectile  with  frontal 
surface  partially  embedded  in  a 
two-layer  target 
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as.  they  relate  to  spheri- 
cal cavity  expansion,  and 
this  quantity  will  he  used 
as  a scaling  function  in 
the  determination  of  the 
axial  force  which  opposes 
the  motion  of  the 
projectile. 


Figure  6.  Projectile  with  frontal 
surface  in  contact  with  both  tar- 
get layers  at  the  same  time 


Denoting  the  pro- 
jectile velocity  by  v and 
the  target  density  by  p , 
the  followirig  ad  hoc  expres- 
sion is  proposed  for  the  ef- 
fective stress  nor- 

mal to  any  portion  of  the 
projectile  frontal  surface 
which  is  in  contact  with 
layer  1: 


a (1)  = 

n 


(1  + 


sin 


ot)a. 


+ vVpi^ 


(21) 


The  subscripts  1 and  2 are  used  to  designate  the  evaluation  of  quan- 
tities in  layers  1 and  2,  respectively.  For  conical  nose  shapes, 
the  quantity  a represents  the  cone  half-angle,  which  is  related  to 
the  nose  length  L and  the  base  diameter  D by 


sin 


-1/2 


(22) 


For  ogival  nose  shapes,  a represents  the  cone  half-angle  at  the 
nose  tip,  which  is  related  to  the  ogive  caliber  radius  (CRH)  and  to 
L and  D by 


sin  o = 


- (ttCRH  - 1)^^^  _ 
2CRH 


Ul/d 


ItL^/D^  + 1 


(23) 


The  evaluation  of  by  equation  20  is  the  same  as  the  evaluation 
of  a for  an  equivalent  spherical  cavity  with  a radius  equal  to  the 
maximiim  radius  of  contact  a^^  between  the  projectile  and  layer  1 
(Figure  5).  Wlien  the  frontal  surface  is  only  partially  embedded  in 
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the  target,  the  equivalent  spherical  cavity  is  centered  at  the  tar- 
get surface,  i.e. 


h = T , z < L 
o * — 


where  T is  the  thickness  of  layer  1 and  z is  the  penetration 
depth.  When  the  frontal  surface  is  fully  embedded  in  the  target, 
the  equivalent  spherical  cavity  is  centered  at  the  base  of  the  nose, 

i.e. 


h^  = T-  z+  L,  L£z£T  + L 


The  evaluation  of  b , k , and  h is  the  same  as  in  Equations  l4- 
l8  with  a = a^  . 


In  accordance  with  Equation  21,  the  effective  stress 
On(2)  normal  to  any  portion  of  the  frontal  surface  which  is  in 
contact  with  layer  2 is 


a^(2)  = (1  + sin  a)Og  + v 


The  evaluation  of  C2  is  the  same  as  the  evaluation  of  a for  an 
equivalent  spherical  cavity  with  a radius  equal  to  the  maximum 
radius  of  contact  a2  between  the  projectile  and  layer  2 (Figure  6), 
This  equivalent  cavity  is  centered  at  the  layer  interface  (h^  = 0), 
and  Equation  20  reduces  to 


= I ^2^ 


1 + In 


Tangential  stresses  are  neglected,  and  the  axial  component 
of  the  effective  normal  stress  is  integrated  over  the  embedded  por- 
tion of  the  frontal  surface.  The  projectile  equation  of  motion  then 
becomes : 


m ^ = -Tra^Cnd)  , z<T 


“ ^ “ “''^(^1  ~ " ^^0^(2)  , T £ z £ T + L 


m ^ = -Tfa^a  (2)  , z ^ T + L 


where  m is  the  projectile  mass  and  t is  time.  The  right-hand 
side  of  the  equation  of  motion  represents  the  toted  eixial  force 
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exerted  on  the  projectile  by  the  target.  . 


Whenever  the  condition  ^ satisfied,  the 

material  in  layer  1 uMgrgoes^cqntained  plastic  flow,  and  Equa- 
tion 21  represents  a"  reasonable  ajjf>roximation  for  Ujj(l)  . However, 
if  this  condition  is  not  satisfied,  then  the  containment  provided 
by  layer  2 may  be  insufficient  to  prevent  a loss  of  cohesion  in 
layer  1,  and  Equation  21  then  represents  an  upper  bound  for  the 
value  of  this  event,  a lower  bound  can  be  obtained  for 

by  introducing  the  following  "plastic  interface  criterion:" 


an(l1 


If  when  b = h (i.e.  when  the  spherical 

plastic  front  located  at  r = b reaches- the  effective  layer  inter- 
face located  at  r = h),  then  Equation  21  is  replaced  by 


0^(1)  = (1  + sin  ct)ag  + vVPgdg' 


(31) 


PENETRATION  AND  PERFO- 
RATION OF  CONCRETE  SLABS 


Canfield  and  Clator  (6)  have  investigated  the  high  speed 
penetration  of  steel  projectiles  into  thick  slabs  of  5000-psi  rein- 
forced concrete,  and  their  experimental  results  are  shown  in  Fig- 
ure 7*  The  projectile  characteristics  are  m = 5.9  kg  , 

D = 76.2  mm  , and  CRH  = 1.5  (ogival  npse  shape,  L/D  = 1.12).  The 
target  properties  are  p = 2.31  gm/cm^  , E = 2k0  kbar  , and 
Y = 0.3^7  kbar  . The  projectile  equation  of  motion  developed  in  the 
present  work  is  integrated  to  obtain  the  theoretical  relation  be- 
tween final  depth  and  impact  velocity  (Figure  7),  based  on  the  as- 
sumption of  a single  semi -infinite  layer  of  homogeneous  concrete. 

The  empirical  results  shown  in  Figure  7 are  obtained  from  a pene- 
tration nomogram  for  concrete  developed  by  the  National  Research 
Council  (NRC)  committee  on  passive  protection  against  bombing  (7). 
The  NRC  nomogram  is  allegedly  accurate  to  within  15  percent,  and  the 
agreement  achieved  among  the  empirical,  experimental,  and  theoreti- 
cal results  is  typicsil  for  cohesive  targets  in  which  the  distance 
from  the  nose  tip  to  the  back  face  of  the  target  is  large  in  com- 
parison with  the  projectile  diameter  (T  - z » D). 


The  results  shown  in  Figure  7 provide  partial  verification 
of  the  applicability  of  the  theory  for  cohesive  targets  composed  of 
a single  semi-infinite  layer.  Now,  in  order  to  make  an  assessment 
of  the  applicability  for  targets  composed  of  distinct  layers,  the 
theory  is  used  to  calculate  the  penetration  and  perforation  of  a 
concrete  slab  of  finite  thickness.  This  represents  an  example  of  a 
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Figure  T*  Penetration  performance  of  76-mm  steel  projectile 
in  a semi-infinite  'slab  of  5000-psi  concrete 

tvo-layer  configuration  in  which  the  first  layer  (concrete)  has  a 
finite  strength  ''.nd  density  while  the  second  layer  (air  or  vacuum) 
has  a negligible  strength  eind  density.  The  concrete  properties  and 
the  projectile  characteristics  are  identical  with  those  used  in  the 
foregoing  calculation  for  a semi -infinite  concrete  slab,  but  in  the 
present  calculation  the  projectile  equation  of  motion  is  integrated 
numerically  to  determine  the  minimum  impact  velocity  required  for 
complete  perforation*  of  a slab  with  a given  finite  thiclmess.  The 
theoretical  results  are  compared  with  empirical  results  in  Figure  8. 
The  empirical  results  are  obtained  from  an  NRC  perforation  nomogram 
for  concrete  slabs  (7)  which  is  allegedly  accurate  to  within 
15  percent. 

Two  sets  of  theoretical  results  are  presented  in  Figure  8, 


* Complete  perforation  means  that  the  projectile  frontal  surface 
passes  completely  through  the  slab,  achieving  a penetration  depth 
z ^ T + L (Figures  5 and  6). 
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Figure  8.  Minimum  thickness  of  5000-psi  concrete  required  to 
prevent  complete  perforation  by  a 76-mm  steel  projectile 


corresponding  to  ceilculations  made  with  and  without  the  plastic  in- 
terface criterion  (Equation  31).  The  theoretical  results  form  a 
fairly  tight  band  about  the  empirical  results,  but  the  calculations 
made  with  the  plastic  interface  criterion  produce  the  best  agreement 
with  the  NRC  nomogram.  The  end  points  of  the  empirical  curve  coin- 
cide with  the  ranges  of  impact  velocity  and  slab  thickness  for  which 
the  nomogram  is  applicable. 


CONCrAJSION 


The  projectile  equation  of  motion  which  has  been  developed 
herein  appears  to  be  applicable  for  cohesive  slabs  and  layered  tar- 
gets in  which  the  thickness  of  a given  slab  or  layer  is  greater  than 
the  projectile  base  diameter  by  at  least  a factor  of  five.  The 
spherical  cavity  expansion  sinalysis,  which  is  used  as  the  basis  of 
the  penetration  theory,  is  inappropriate  when  the  diameter  of  the 
projectile  is  comparable  wit*  the  dimensions  of  the  .target.  The 
analysis  of  the  penetration  of  thin  target  layers  (T  ^ D)  requires 
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a different  conceptual  approach  whicTi  lies  beyond  the  scope  of  the 
present  work.  The"  theory  can  be  used  to  calculate  the  instantaneous 
deceleration,  velocity,  and  position  of  a rigid  cylindrical  projec- 
tile (with  conical  or  ogival  nose  3hape)  when  the  direction  of 
motion  is  coincident 'with  the  pirojectile  axis  of  symmetry  and  normal 
to  the  target  surface.  The  de+a-.ls  of  the  target  motion  are  not 
predicted  by  the  theory,  which  relies  on  a highly  idealized  model  or 
the  projectile/target  interaction.  Two-dimensional  finite- 
difference  solutions  represer.t  the  o,  ly  means  pi  -sently  available 
for  assessing  the  details  of  the  target  bch'^tor  during  the  pene- 
tration process. 
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